Fermi- surface induced modulation in an optimally doped YBCO superconductor 
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We have observed a Fermi- surface (FS) induced lattice modulation in a YBCO superconductor 
with a wavevector along CuO chains, i.e. q^ = (0,(5,0). The value of ^ ^ 0.21 is twice the Fermi 
wavevector (2kF) along b* connecting nearly nested FS 'ridges'. The qi modulation exists only 
within 0-vacancy-ordered islands (characterized by qo=(|, 0, O)) and persists well above and below 
Tc. Our results are consistent with the presence of a FS-induced charge- density wave. 

PACS numbers: 74.72.Bk, 61.05.cp, 71.45.Lr 



YBa2Cu307_aj (YBCO) compounds possess a unique 
structure where one-dimensional CuO chains are inter- 
leaved with two-dimensional Cu02 bi-layers, from which 
superconductivity appears to emerge. These chains not 
only act as the charge reservoirs for the Cu02 planes, 
but also provide the opportunity for the interplay of elec- 
tronic properties of systems with different dimensionality 
in close proximity to each other. Various studies [1] |2l |3] 
have revealed that the chains exert profound influence 
on the transport properties in both normal and super- 
conducting states. Local density approximation (LD A) 
calculations [4] have indicated sections of Fermi surface 
(FS) that are open electron "ridges" (|| FX) orthog- 
onal to b* {i.e. normal to the CuO chain direction, 
see Fig. [l]). These rid ges are due to CuO-chain derived 
electronic states, and have been confirmed by positron 
annihilation [5^ and photoemission[6 studies. Intrigu- 
ingly, near the zone center these ridges remain as parallel 
sheets [7 along c* forming a near nesting condition with 
an incommensurate wavevector qic = (0,2/ci? ~ 0.22,0) 
(in reduced lattice units; k^? is the Fermi wavevector), 
suggesting that the CuO chains may be on the verge of an 
electronic instability {e.g. a charge-density wave, CDW). 

Although local probes such as NMR/NQR 
measurements [8 have reported on the formation of 
a CDW on the CuO chains, the results of different 
groups have been at odds with each other. Observations 
of electronic modulations of the chains using STM[TOl[TT] 
have seemed promising because the observed periodicity 
was consistent with the expected qic for an optimally 
doped YBCO, but doubts have been raised [12] on the 
CDW interpretation due to the dispersive behavior of 
the modulation wavevector on energy. Furthermore, 
it remains unclear whether a static CDW on the ID 
CuO chains as was suggested by STM workfOT, TT 
forms in the bulk or not. Inelastic neutron scattering 
(INS) measurements [13 have also been used to infer 
the presence of charge fluctuations with wavevector of 



(0, ~ 0.23,0). These measurements, however, could 
not distinguish between a and b directions due to the 
twinned samples used in the experiment, adding an 
additional directional ambiguity. 
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Figure 1: Left: K scans (T~7K) showing shoulders on both 
sides of the central qo satellite peaks (lines are from Gaussian 
profile as a guide to the eye). Vertical dashes indicate the 
shoulders displaced along K by ±0.21 from their respective 
qo satellite. Right: A schematic showing near nesting of the 
'ridges' (purple) with nesting vector along b* (red) O [6l [7] . 

In this Letter, we present x-ray scattering observa- 
tions of lattice modulations with a wavevector qi = 
(0, ~ 0.21, 0) along the CuO chain direction. Atomic dis- 
placements associated with them are parallel to the or- 
thorhombic a axis. They co-exist with the 0-ordered 
ORTHO-IV phase [Ml [15] islands characterized by qo = 
(±|,0,0), rather than forming homogeneously through- 
out the crystal. In reciprocal-space scans along b* (K 
scans. Fig. [T]), the qi super-lattice peaks manifest them- 
selves as shoulders on either side of the central qo satel- 
lites. Below we summarize our previous finding, then 
show the coupled nature of these two modulations. A 
quantitative modeling of the x-ray diffuse scattering pat- 
tern in the H-K plane is used to establish the source of 
the shoulders. 

A high-quality detwinned crystal of optimally doped 
YBCO7-X {x = 0.08 and = 91. 5K) was used in our 
study, which eliminates any directional ambiguity. We 
note that the presence of qo and qi modulations was 
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verified in a twinned crystal as well which will be pub- 
lished elsewhere. High energy (E = 36keV) x-ray dif- 
fuse scattering experiments were performed on the X- 
Ray Operations and Research 4-ID-D beam line at the 
Advanced Photon Source (APS). The choice of high x- 
ray energy makes this a truly &t^/A;-sensitive study [M]. 
At this composition, first principles band calculations [15] 
have shown that the oxygen vacancies, located in the 
chains, tend to cluster together to form a minority phase 
with a 4-unit-cell superstructure, corresponding to a se- 
quence of three 0-full and one 0-vacant chains (ORTHO- 
IV phase denoted by < 1110 >, see Fig. |4|, respec- 
tively, the formation of which was confirmed in a pre- 
vious X-ray scattering study [14 . That study revealed 
that, in the ORTHO-IV phase islands, coupled displace- 
ments of Cu, Ba, and O are induced, respectively, in 
the a-c plane with anisotropic correlation lengths along 
each of the three crystallographic axes, producing broad 
satellite peaks centered at qo = (±^,0,0) on both sides 
of Bragg points. Embedded in the average orthorhombic 
crystal, the ORTHO-IV phase islands induce long-range 
strain in the lattice which gives rise to the Huang diffuse 
scattering (HDS)^ around Bragg peaks (Fig. |2|a)). 

If lattice modulations associated with CDWs with 
propagating vector qic are present throughout the whole 
crystal, they should be observed by X-ray scattering as 
peaks (or rods) centered at (m, n, L) ± qjc {m and n 
are integers) in the H-K plane of the reciprocal space, 
depending on their inter-chain correlations. A thorough 
investigation in the vicinity of (m, n, L) ± qic for dif- 
ferent integer m and n values showed no indication of 
such modulations along b* direction as claimed by the 
INS measurements [13]. Instead, clear shoulders were ob- 
served at qo ±qi, qi = (0, (5, 0) with S ~ 0.2, as shown in 
Fig. [l] It is interesting to notice that K-scans through 
different qo peaks show a definite correlation between the 
intensities of the shoulder peaks with the center ORTHO- 
IV peaks, both having the same correlation lengths and 
obeying the same parity relations [14 . The parity rela- 
tions impose stringent phase correlations between atomic 
displacements such as those in the ORTHO-IV model [14 . 

Furthermore, careful measurements of the ratio of the 
satellite and the shoulder intensities as a function of in- 
creasing temperature found it to be constant in the range 
of 7-300K. As pointed out in Ref. [14], the temperature 
dependence of the qo peaks distinguishes the ORTHO- 
IV phase islands from the average lattice which has less 
rapidly varying Debye- Waller factors. The fact that the 
intensities of the shoulders obey the same temperature 
dependence as the center qo peak strongly indicates that 
they are intimately related. Therefore, the qi modula- 
tion must coexist with ORTHO-IV islands rather than 
forming homogeneously throughout the whole crystal. 

Although a commensurate 0)-type super- 

structure in YBCO had been reported in the early 
literature [17 , it was ruled out to be due to O- vacancy 



ordering in later studies [T8^. Such superstuctures should 
give rise to a quartet of (±|,±|,0) satellites of Bragg 
peaks, which is not observed here, suggesting that they 
were peculiar to the samples used in the electron diffrac- 
tion work, and are not a general property of the bulk. 
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Figure 2: (a): Colormap of the measured diffuse scattering 
around (4, 0, 0) and (5, 0, 0) in H-K plane, (b) Model calcula- 
tions of the 2D scattering pattern. Asymmetric HDS around 
(4, 0, 0), qo satellite, and qi shoulders are all reproduced. 
The dashed lines indicate the line cuts presented in Fig. [3] 

In order to identify the origin of the shoulders as shown 
in Fig. [l] two-dimensional scans around several Bragg 
peaks were performed and the diffuse scattering inten- 
sity pattern around (4, 0, 0) and (5, 0, 0) is shown in Fig. 
|2|a). The salient feature of this pattern is the profound 
difference in the symmetry of scattering around the two 
Bragg peaks. While the dominating feature around the 
(4, 0, 0) Bragg peak is a strong nearly fourfold 'bow-tie'- 
shape HDS pattern albeit asymmetric, it is entirely miss- 
ing around (5, 0, 0). Rather, there is a strong 'elliptical' 
pattern centered at the (5.25, 0, 0) ORTHO-IV satellite 
with additional scattering bulging out along the K axis, 
which appears as shoulders in K scan through (5.25, 0, 
0). The intensity near (3.50, 0, 0) is from the tail of the 
strong qo-satellite at (3.25, 0, 0). In particular, the inten- 
sity asymmetry about the K axis through (4, 0, 0) is due 
to the interference scattering (IFS) from the ORTHO-IV 
phase islands and the lattice-strain fields, respectively. 
Such interference effects have been known in the study 
of disorder and extensively discussed in Refs.[19j. 

The presence of both HDS and IFS has to be taken into 
account in discussing the origin of the shoulders. The 
diffuse scattering pattern around (4, 0, 0) shown in Fig. 
|2|a) suggests that the shoulders in the K-scan through 
(4.25, 0, 0) in Fig.[3];a) might be due to the lobes of HDS, 
extending roughly along [110] directions in H-K plane. 
But this scenario is unlikely to account for the strong 
shoulders of the (5.25, 0, 0) satellite given the absence of 
similar HDS around (5, 0, 0) Bragg peak. 

The absence of HDS around (5, 0, 0) can be understood 
if one examines the general nature of HDS [19 . As a long- 
range elastic distortion of the lattice, the strain fields 
give rise to HDS in the vicinity of a given Bragg point in 
proportion to the square of the unit cell structure factor 
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modulus, |F(G)|^ = lEfe Ae"^'^^^e^^■^M^ where G 
is the Bragg point[19 . For YBCO6.92, l^fHyT - 0.02. 
With this ratio, we estimate the HDS at (5.25, 0.2, 0) to 
be about ~ 20 counts above background for the K-scan 
in Fig. [sj assuming the shoulder at (4.25, 0.2, 0) is from 
HDS. The tail of the qo peak is estimated to give ~ 20 
counts at (5.25, 0.2, 0) from Gaussian fitting. As a result, 
the interference term, IFS, can be at most ~ 40 counts 
at (5.25, 0.2, 0). Both the HDS and IFS are negligible 
compared with the observed strong (5.25, 0) shoulders 
which are about ~ 2000 counts. Thus, neither HDS nor 
IFS can account for the shoulders. They are the finger 
print of a lattice modulation along the b axis. 

We first discuss the X-ray diffuse scattering without 
considering the b-direction modulation. With the pres- 
ence of the ORTHO-IV phase islands and their surround- 
ing strained lattice, the total diffuse scattering intensity, 
/diffuse (Q) 5 can be written as in Eqn. (1). A linear ap- 
proximation is made assuming that the lattice distortions 
due to the strain are small. 



^diffuse(Q) < I fk{Q)e 



jl'k 



h^^Q-v,ifeA(Q)e-^^(^)e-^Q-^-p> (1) 



Ik 

j is the index for the 4a x lb x Ic supercell[14 of an 
ORTHO-IV minority phase island, ^ji' k lattice 
distortion vector of the k-th atom of the / -th unit cell 
inside the j-th supercell. Outside of the island, u^-^'/^ is 
zero and the first term vanishes, v/^ is the lattice dis- 
tortion of the k-th atom of the l-th unit cell due to the 
strain and the summation of the second term is over the 
whole crystal. fk{Q) and e"^''^^^ are the atomic scat- 
tering factors and Debye- Waller factors (both dynamic 
and static), respectively. The modulus squared of the 
first term gives rise to the ORTHO-IV phase satellites, 
the modulus squared of the second term describes the 
HDS from lattice strain and the interference of these two 
terms gives rise to IFS. The additional thermal diffuse 
scattering was calculated using the shell model by fit- 
ting to the measured phonon dispersion curves [20] and 
subtracted from data except from those in Fig. [l] 

The scattering from different contributions dominate 
different regions in reciprocal space. For example, the 
HDS has a minimum at all qo positions and gives neg- 
ligible contributions around Bragg peaks with odd par- 
ity. This allows us to decouple most of the parameters 
and fit them separately. In particular, in Ref.[14 , u^^'/^ 
has been fitted to the integrated intensities of a set of 
G ± qo peaks alone (G is a reciprocal lattice vector) 
which were extracted by a careful line fitting procedure, 
and an ORTHO-IV displacement model was constructed. 

Now, we need to incorporate the b-direction modula- 
tion into Eqn.(l). As discussed above, the new modula- 
tion does not form homogeneously throughout the whole 
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Figure 3: Solid lines are the fitting results. 6 and d are de- 
termined to be 0.21 and 0.9, respectively. The K scan in (a) 
across (4.25, 0, 0) is over estimated by ^ 13% mainly due to 
the IFS contribution and the lower side of the H scan in (b) 
is pushed up by the tail of strong (3.25, 0, 0) satellite peak. 



crystal. Rather, a qi = (0, ^, 0)-type wavevector acts as a 
b-direction periodic modulation of the displacements of 
all the atoms involved in the ORTHO-IV phase [M]. We 
incorporate this new modulation with a simple extension 
of the ORTHO-IV displacement pattern by replacing u^s 
in Eqn. (1) with t's as follows: 



Uk = u//e[l + dcos (qi -R, 



Ik 



(2) 



where / is the index of the unit cells inside the ORTHO- 
IV phase island and k is the atom index in that unit 
cell. Since the u//., as reported in Ref.[14], lie in the a-c 
plane, the second term of Eqn. (2) represents an addi- 
tional transverse modulation along the b-direction. The 
second term contributes to the qi shoulders and has small 
effects on qo peaks. Thus fitted values of uik remain close 
to those reported in Ref. [14 within the errors introduced 
by the IFS term included in this work. 

To determine the values of 6 and the experimen- 
tal data is fitted with the ORTHO-IV model, together 
with the new lattice modulation as in Eqn. (2). With 
this model the 2D pattern can be reproduced remark- 
ably well as shown in Fig. |2|b). Fig. [s] gives a better 
quantitative view of the data fit along important lines 
indicated in Fig. |2j Thus, our experimental observation 
and quantitative analyses establish the existence of b- 
direction lattice modulation with qi = (0, ^ ~ 0.21,0) 
within the ORTHO-IV phase islands in optimally doped 
YBCO. qi is very close to qic, strongly indicating that 
this modulation is induced by the FS. 

The idea that FS nesting can infiuence short-range or- 
der (SRO) and manifest itself in the X-ray diffuse scatter- 
ing was discussed [21] in the study of order-disorder tran- 
sitions in binary alloys such as Cu-Au and was used sub- 
sequently in investigations [22] of SRO correlated micro- 
domain structures within a disordered matrix. In the 
context of the ORTHO-IV phase any one of the four per- 
mutations of < 1110 > can form as anti-phase (AP) do- 
mains with periodicity along the b axis determined by qi 
giving rise to symmetric shoulders of qo, along K axis. 
However, structure-factor calculations based on various 
combinations {e.g. two < 0111 > and three < 1110 >; as 
illustrated in Fig.|4| assuming no amplitude modulations 
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Figure 4: Left: ORTHO-IV displacement pattern (a6-plane) 
for the chains [14]. qi further modulates amplitude of these 
displacements sinusoidally. Right: One possible AP domain 
with no modulations of displacement amplitudes. Small cir- 
cles (red): Cu; Big circles (blue): Oxygen; Squares: Vacancy; 
Dashed line: Domain boundary; Shaded area: 4X5 super-cell. 



of the ORTHO-IV displacement pattern [14 did not yield 
relative intensities of qo and qi satellites consistent with 
the data. Although we can not rule out a more subtle 
arrangements of AP domains, a simple transverse modu- 
lation of a mono-domain ORTHO-IV pattern, as modeled 
in Eqn. (2), adequately accounts for the data. 

It is important to note that the shoulders get con- 
tributions from all the transversely displaced atoms in 
the ORTHO-IV islands (predominantly the Ba atoms), 
rather than from the chains alone. In our view, nested 
'ridges' (or 'sheets') of the FS simply induce a qi ampli- 
tude modulation of the ORTHO-IV displacement pat- 
tern with a concomitant presence of a CDW on the 
CuO chains. Since the ORTHO-IV displacement pat- 
tern is a static relaxation of the lattice [15], we expect 
the qi modulation to be static as well. The formation 
of qi modulations only inside the ORTHO-IV phase is- 
lands rather than throughout the average lattice seems to 
disagree with surface-sensitive STM experiments [TOl [TT] 
which have reported well-patterned electronic corruga- 
tions in regions away from oxygen vacancies. This dis- 
crepancy may be due to energetics of CDW formation 
on chain-terminated surfaces as opposed to that in the 
bulk. The LDA calculations [4 indicate that the CuO 
chains in the bulk are on the verge of a CDW instability 
due to FS nesting, but whether such a CDW can form 
or not depends on energy gain of carriers relative to the 
energy costs due to the lattice deformations. The lack of 
any peaks (or rods) in the vicinity of (m, n, L) ± qi in 
present study indicates that the average lattice may be 
too rigid for a CDW to form. Instead, it seems that a 
static CDW can only be established in favorable regions, 
such as the cleaved surface with CuO chains to be the 
first layer, as studied with STM, or in other elastically 
"softer" areas. A lattice softening introduced by oxygen 
vacancies is suggested by STM observations [lOl HTJ [12] 
of enhanced CDW amplitudes in their vicinity. X-ray 



studies [11123] have also shown that ORTHO-IV islands 
are softer than the average lattice with structurally co- 
herent chains, which in principle can favor the CDW for- 
mation. Since the displacements are transverse to qi and 
are three dimensional in nature, a straightforward inter- 
pretation based on a canonical longitudinal CDW seems 
precluded. However, we can speculate that ORTHO-IV 
islands may act as 'pinning' templates [23j for dynamic 
charge fluctuations [24], calling for further studies of FS 
effects in YBCO. 

In conclusion, in optimally doped YBCO, an incom- 
mensurate transverse lattice modulation with wavevec- 
tor qi = (0, - 0.21, 0) is observed inside the ORTHO-IV 
phase islands. This modulation is induced by FS nesting 
associated with ID CuO chain subsystems and forms in- 
homogeneously in the lattice. Since the Cu02 planes are 
also involved in the qo ± qi modulation, FS effects must 
affect the planes as well. 
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